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Abstract 

We study homogenization problem for the stationary Maxwell sys- 
tem. It is supposed that the magnetic permeability and the dielectric 
permittivity locally close to fast-oscillating (with respect to some small 
parameter) periodic functions which can change the form on rather big 
distances. An asymptotic behavior to solutions of the Maxwell sys- 
tem outside of its spectrum is obtained. We also describe asymptotic 
behavior of resolvent with control of the remainder in terms of some 
appropriate operator norms. 

1 Introduction 

The work is devoted to study the homogenization problem for the Maxwell 
operator. Roughly speaking it is supposed that the magnetic permeability 
and the dielectric permittivity locally close to fast-oscillating periodic func- 
tions which can change the form on rather big distances. Problems of this 
kind for operators with purely periodic coefficients were intensively studied 
by many authors [U, El EJ EJ [71 El The Maxwell operator with purely pe- 
riodic coefficients was in details studied in a series of works of M.Sh.Birman 
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and T.A.Suslina, last results can be found in work [U]. However, homog- 
enization problems concerned to operators whose coefficients aren't purely 
periodic are investigated less in details [H [5] . Moreover, it seems like, the 
Maxwell operator with coefficients locally close to fast-oscillating periodic 
functions remains almost not studied. 

Let us denote by L2(M^; C^) the L2 class of C^- valued functions in and 
by H^{M?\ the corresponding Sobolev classes of order p, where p G N. 
Let us put 

J ={u:ue L2(M^; C^), div(M) = 0}. 
Here equality div(n) = is understood in sense of distribution theory 

diY{u) = 0^ J{u,Vw)dx = V weH\W^;C^), 

where (■, ■) is the standard scalar product in C^. 

Let us consider the model described by the Maxwell operator 



u\ _ f zrot(yU ^■)\ I'u 



where a and fi are some matrix-valued functions which will be described in 
detail lately. The Maxwell operator 2Jl acts in the space J© J on the domain 

j:)omm= {u,v:u,v e H\R^,C^), div(M) = 0, div(w) = 0}. 

Now let us describe a and /i. Let F be a lattice in and Q be the 
elementary cell of the lattice F. Let coefficients a and fi have the form 

{ x\ ( X 

a = a \ x, — j , fi = fi yx, — 

where e is a small positive parameter. Suppose that the matrix-valued func- 
tions a and fi satisfy the following assumption. 

Assumption 1.1 

1) a and fi are real and positive- definite (3 x 3) matrixes; 

2) a, fie C~(M^M3); 

3) a{x,y) and fi{x,y) are T-periodic with respect to y; 

4) a{x,y) = I and fi{x,y) = I for any x ^ Br = {|a;| < R} and y G M^, 
where I is the identity matrix and R is a positive constant. 
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Denote by £Dt(e) the Maxwell operator with coefficients a and [i satisfying 
assumption 11.11 The Maxwell operator 9Jt(£:) is closed in J © J with the 
standard scalar product and selfadjoint in J © J with scalar product of the 
form 

We study an asymptotic behaviour as e — )■ of the resolvent (^{e)—E)^^ 
for Imi? > 0. We also investigate asymptotic behaviour of solutions to the 
Maxwell system 

where \va.E > and 

r, r e c°°(M^ c=^) n L2(M^ c^), div(/") = o, div(/^) = o. (2) 

The main results of the work are contained in theorem 13.11 and corollary 



2 Formal asymptotic solutions 

This section is devoted to the constructing of formal asymptotic solutions, 
as e — )■ 0, of the system ([T]). To separate the slow and fast dependencies on 
the argument we seek the solution of system ([T]) in the form 



U{x,e)\ _ /'u{x,x/e,e) 
V{x,e)J \y{x,x/e,e) 

It is easy to see that the following lemma holds. 

Lemma 2.1 Let functions u and v satisfy equations 

uix,y,e)\ _ ( f'ix) 



(1) 



e ^ divyu{x,y,e) = div^u{x,y,e), e ^ divyv{x,y,e) = div^v{x,y,e), (3) 
where 



u{x,y)\ ^ / irot^ifi ^{x,y)v{x,y)) 
v{x,y)) \-iiot^{a-^{x,y)v{x,y)) ) ' 



fu{x,y)\ ^ / iiotyifi ^{x,y)v{x,y)) 
^ \v{x, y) J \-i roty{a-^{x, y)v{x, y))^ 

and the operators mt^, dhr^ and loty, div^ act upon a variables x and y 
respectively. Then a function 

U{x,e)\ ^ ^u{x,^,e) 
V{x,e)J \v{x,j,e) 

satisfies equation ([1]) and conditions dw{U{x,6)) = 0, dw{V{x,e)) = 0. 

Proof. It can be easily proved by using direct calculations. □ 

Thus, the constructing of formal solutions of system ([T]) is reduced to the 

constructing of formal solutions of equations ([2]), ([3]) 

The formal solutions of equations (|2]), ([3]) are seeked in the form 

uix, y, 6)\ _sr^ n f Unix, y)\ , . 

2^^ vJx^v) ' 



v{x,y,e)J ^ \Vn{x,y) 



n>0 



where functions Un{x,y) and Vn{x,y) are F-periodic with respect to y. 

Lemma 2.2 Suppose Un, satisfy the following recurrence system of equa- 
tions 

<l\Vy{un{x,y)) = -div^{un-iix,y)), divy{vnix,y)) = - div^(f„_i(x, y)), 

(6) 

where m_i = and f_i = 0. Then (jl]) formally satisfies problem ([2]), ([3D. 

Proof. Substituting expansion (jl]) in system ([2]), ([2]), and comparing the 
coefficients corresponding to the equal powers of e, we obtain recurrence 
system ([5]), ([6D. □ 

For detailed description of solutions to system ([5]), ([6]) we recall the fol- 
lowing well known result. 

Lemma 2.3 There exist 3 linearly independent T -periodic (with respect to 
y) solutions to the system 

miy{C{x,y)) = 0, 
diVy{a{x,y)Cix,y)) = 0. 



(k{x,y) = ek + Vyfk{x,y), k = 1,2,3, 
ei = (1,0,0)*, 62 = (0,1,0)*, 63 = (0,0,1)*, 
where ipk{x,y) is any T -periodic (with respect to y) solutions to the equation 

divy{a{x, y)Vyipk{x, y) + a{x, y)ek) = 0. 

There exist 3 linearly independent V -periodic (with respect to y) solutions 
to the system 

ioty{i{x,y)) = 0, 
divy{fi{x,y)^{x,y)) = 0. 

^k{x,y) = Ck + VyiJk{x,y), k = 1,2,3, 

where ilJk{x,y) is any T -periodic (with respect to y) solutions to the equation 

diVy{n{x,y)Vyipkix,y) + fi{x,y)ek) = 0. 

Moreover, Ck,^k e C°°(M.^,C^) and Ck{x,y) = e^, ik{x,y) = Ck for x ^ 
Br, where k = 1,2,3. 

Proof. We only need to prove that (k{x,y) = Ck, C,k{x,y) = Ck for x ^ Br, 
where k = 1,2,3. It is easy to see that a{x,y) = fi{x,y) = / for x ^ Br. 
Hence ipk, ipk satisfy the following equations 

Ayipkix,y) =0, Ay^kix,y)=0, x ^ Br, k = 1,2,3. 

Since ipk{x,y), ipk{x,y) are F-periodic functions (with respect to y), we see 
that ipk{x,y), ilJk{x,y) are constant for x ^ Br. This implies the necessary 
assertion. □ 

Consider recurrence system (|5]), ([6]) for n = 

"^y = 0' diYy{u,{x,y)) = diYy{v,{x,y)) = 0. (7) 



y \vQ{x,y) 

Lemma 12.31 implies that general solution of ([7]) has the following form 

3 3 

uq{x, y) = a{x, y) ^ al{x)Ck{x, y), vq{x, y) = /i(x, y) ^ bl{x)^k{x, y), 

k=l k=l 

where a°(x) and b^{x) are arbitrary smooth functions. 
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Recurrence system for n = 1 has the from 



r{x) 



divy{ui{x,y)) = -div^{uo{x,y)), divy{vi{x,y)) = -div^{vo{a 
Let us rewrite it in the following from 

divj^(fi) = -div^{vo), 

rotj^(a;^-'^Mi) = if" + iEvq — rot3;(a^"'^Mo), 
diVj^(Mi) = -div^(Mo). 
Suppose that the following conditions hold 

/ {-if - iEuq - rot,.(^-it;o), a) dy = Q, = 1, 2, 3, 
n 

J div^(t;o) dy = 0, 
n 

then problem (j9]) can be solved. It is easy to check that 

1^ J i^iix, y),^2ix, y),^3ix, y)y dy = I, 



'{nx),Ux,y))\ 

{nx),Ux,y))^ 

{nx),Ux,y))J 



\n\nx). 



Similarly, 



{uo{x, y),ik{x, y))y = \n\ J] A^p(x)aJ(x), 



( {uo{x,y),^i{x,y))\ 

{uo{x,y),^2{x,y))y 
\{uo{x,y),^3{x,y))yJ 



|fi|A"(x)a°(x), a'\x) 



0'2ix) 
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where 



= -^^{(^{x,y)Cp{x,y),Ck{x,y))^. (12) 

From f[T2]) it follows that A"(x) is real and positive-definite (3 x 3) matrix. 
Using direct calculations, we get 

3 

{Tot^{fi-^{x,y)vo{x,y)),^k{x,y))y = ^ (rot^(6°(x)^p(a;,y)),^fc(x,?/))^ = 

p=i 



(rot^(^p(x,?/)),^fc(x,?/))^ 



p=i 



p=i 



It is easy to see that 



{Tot^{^p{x,y),^k{x,y))y = 0, 
[^pix,y) X ^kix,y)] dy = \n\ [cp X Cfc] . 



Therefore, 



{Tot^vo{x,y),^k{x,y))y = \n\Y{ [cp X Cfc] , V^(6°(x))) 

p=i 



I (rot^ t^o(a;,l/),6(x,?/)) \ 

(rot^ VQ{x,y),i2{x,y))y 
y(rot^ Vo{x,y),^3{x,y))^J 



'hl{x) 

\n\mi^lP{x), b%x) = I 6° (a;) 



Thus, the first three equations from (|TT1) can be rewritten in the following 
form 

2rot^6°(a;) - EA"(a;)a°(x) = /"(x). 
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The last equation from fllip can be reduced to the form 

div,(A^(x)6°(x))=0, 

where 

(13) 

From f|T3|) it follows that A^(a;) is real and positive-definite (3 x 3) matrix. 

Similarly, one can check that solvability conditions for problem ( 1T0|) can 
be written in the form 

-iiot^a'^ix) - EA''{x)b\x) = r{x), 

div^.(A"(a;)a°(a:)) = 0. 
Let us introduce notations 

uo{x) = A"(x)a°(x), vo{x) = A''{x)b%x). 

Finally, the solvability conditions for ([8]) have the following form 

diVa;(Mo) = 0, div^(i;o) = 0, 

(14) 

General solution of (IHl) can be represented in the form 

3 

ui{x,y) = ui{x,y) + a{x,y)'^al{x)Ck{x,y), 

k=l 
3 

vi{x, y) = vi{x, y) + ^{x, y) hl{x)ik{x, y), 

k=l 

where ui{x,y), vi{x,y) are partial solutions of ([8]) that are orthogonal to 
solutions of homogeneous equations. 

Let us check that Ui{x, y) = and Vi{x, y) = for x ^ Bji. It is easy to 
see that 

a{x, y) = fiix, y) = A"(2;) = A''(a;) = /, ^0(2;) = a°(a;), vo{x) = 6°(x), 
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for X ^ B^. Taking into account (fT^ . we see that for x ^ -B/j the following 
equations hold 




(15) 



TOty{Vl) = 0, 

divy{vi) = 0, 

Recalling that ui{x,y), Vi{x,y) are orthogonal to solutions of homogeneous 
equations, we get that Ui{x, y) = and Vi{x, y) = for x ^ Br. 

Now let us consider recurrence system ([S]), (jH]) for 2. General solution 
of this system can be represented in the form 



Un{x, y) = u„{x, y) + y) al{x)Ck{x, y), 



k=l 



(16) 



Vn{x, y) = f„(x, y) + /i(x, y) J] bl{x)^k{x, y), 



k=l 



where Un{x,y), Vn{x,y) are partial solutions of ([S]), ([S]) that are orthogonal 
to solutions of homogeneous system. Coefficients a^{x) and b'^{x) can be 
defined from the solvability conditions of system ([5]), ([6]) with n substituted 
by n + 1, which can be written in the form 



where 



(m-E) 



divjur. 



Vn 



Gl{x) 



diVxfWr, 



Gl{x), 



Un{x) = A"(a;) 




Vn[X 



A^ix) 




F^ix') = 



/{f:{x,y),Ux,y))y\ 

{fnix,y),^2ix,y))^ , 
\{fni^,y),^3{x,y))J 

fni^^y) = Evnix,y) + iTot^{fi''\x,y)vnix,y)), 
/{rnix,y)Xiix,y))\ 

{fnix,y)X2{x,y))y , 
\{fn{x,y)X3{x,y))J 

f 'ni.x,y) = Eurr{x,y) + z rot^. (a~^ (x, |/)Ct„(x, ?/) ) , 



1 

W\ 



(17) 
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Glix) = j div^{unix,y))dy, Gl{x) = j div^{vn{x,y)) dy. (18) 
Q n 

Let us collect the obtained result in the following theorem. 

Theorem 2.4 Suppose lm{E) > 0, e > 0, assumption 11.11 and conditions 
(I2l) hold. Then Maxwell system ([T]) has asymptotic solution of the form 



U{x,e] 
V{x,e] 



Ur}. [ X 



Vr,. [X 



(19) 



where Un and f„ can he represented in form f|T6|) . The corresponding coeffi- 
cients can he founded from recurrence equations ([5]), ([6]), and f|T7|) . 
Moreover, Un{x, y) = and Vn{x, y) = for x ^ and n ^ 0. 

The following lemma describe important estimates for the components of 
formal solution ffT9l). 



(20) 

(21) 

(22) 
(23) 



Lemma 2.5 Suppose that the statements of theorem^^ hold. Then 

d\k\ 



max max 



max max 

0!£|fc|^2 yen 



Qyk 

d\k\ 
dy 



Un{-,y,e) 



n—s—1 



<C\\f\\Hs 



for ^ n ^ s + 1, where the constant C does not depend on f 



_n— s— 1 



Proof. The statement of the lemma for n = follows from flMj) . Using 
recurrent system of equations ([S]), ([0]), and f[T7|) . the statement of the lemma 
for n ^ 1 can be easily proved by induction. □ 
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3 Asymptotic expansion of the resolvent 



In this section we describe asymptotic behaviour of the resolvent — 
E)^^ for a small e. To do that we consider a partial sum for formal series 
ffim of the form 



Un{x, x/e, e) 
VN{x,x/€,e) 



N 
n=0 



Vn{x,x/e)J ' ^' 



where N ^ 1 and 6 = (5", 5^)* is an auxiliary function. Here we also suppose 
that vn{x) = un{x) = 0. This implies that UN{x,x/e) = VN{x,x/e) = for 

X ^ Br. 

It easy to see that 



m^) - E) 



N 



Vn 

N 



r{x) 
rix) 



+ 



- E) 



+ e''{mie)-E)5{x,e), (1) 



y=x/e 



VNix,y) 

div (Un) = £^ [div,.(M;v(a;, y))]y=,/, + e"" div(5"(x, e)), (2) 
div {Vn) = £^ [div.(^^^(x,y))],=,/, + £^div(5^(x,£)). (3) 

Now we chose 5 such that {UniVnY € DomOJl. To do that we find an 
appropriate solution to the following equations div(f//v) = 0, div (V)v) = 0. 
First equation can be rewritten in the form 

div(5"(x, e)) = g{x, e) = - [div^(Miv(a:, y))]y=^/e ■ (4) 

Let us consider the following solution of equation (j4]) 

1 

5'^{x,e)=x J g{tx,e)t'^ dt. 



Recalling that supp g C Br, we obtain 

l'^"(-,^)lllfi«3,c3) ^ (i?' + 6) / / \g{tx,e)\t^dt] dx+ 




+ 2R^Y^ J I J \gj{tx,e)\t^dt\ dx, 
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where gj{x,e) = dx.g{x,e). Using simple calculations, we get 



\g{tx,e)\t'dt\ dx^\\g{;e 

R3 \0 

1 V 2 

3 



\gj{tx,e)\rdt\ dx^Wgji- 



e 



>C))' 



k = 1,2,3, 



and hence 



\\n,e)\\ 

+ 2R'J2M-^^)\\liR^c)^{^R' + m9{-,e)\\% 
Using lemma 12.51 and definition of g we obtain 

In the same way one can define 6^ such that div (Vat) = and 
||5"(-,^)bi(R3,c3) ^ {2R' + 6)s~'\\f\\HN+i 



(5) 
(6) 



Since {uiy{x,y),VN{x,y)y G DomOJt, we can apply the resolvent 
E)-^ to the both sides of © for lm{E) > 



where 

h{x, x/e, e) 



- E) 



UN{x,y) 



+ {m{e)-E)6{x,e). (7) 



y=x/e 



Using lemma 12. 5[ ([5]) , ([6]) , and ([7]) we get 



me) - EY 



1 fnx) 
r{x) 



N-2 
n=0 



Un{x,x/e) 
Vn{x,x/e) 



L2(IR3,C6) 

r{x) 
r{x) 
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for N ^ 2, where the constant Ci does not depend on /"(x) and f'"{x). 
Finally, we proved the following theorem. 

Theorem 3.1 Suppose \m{E) > 0, e > 0, assumption 11.11 and conditions 
hold. Then estimate (|H]) holds for N ^ 2, where the constant Ci does not 
depend on /"(x) and f"{x). 

Corollary 3.2 Suppose lm{E) > 0, £ > 0, assumption 11.11 and conditions 
(H hold. Then 



{m{e) - E)-^ - e(£)(aji - e)-^- h%m.\ c^) ^ L2(M^ c^^ 



where constant C\ does not depend on e and Q{e) is a multiplication operator 
on the matrix-valued function 

a(x,x/£)Z"(x,x/e)A"(x)-i \ 

/i(x,x/e)Z^(x,x/£)A^(x)-7 ' 

^" = (Cl,C2,C3), ^^ = (^1,6,6). 

Proof. For the proof it is sufficient to apply theorem 13. II for = 2. □ 
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